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We show in this letter that estimates of the local emissivity of 7-rays in the GeV-TeV range 
suffer uncertainties which are of the same order of magnitude as the current Fermi results. Primary 
cosmic-ray fluxes, cosmic-ray propagation, interstellar helium abundance and 7-ray production cross- 
sections all affect the estimate of this quantity We also show that the so-called nuclear enhancement 
factor - though widely used so far to model the 7-ray emissivity - is no longer a relevant quantity 
given the latest measurements of the primary cosmic ray proton and helium spectra. 
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INTRODUCTION 

The study of the Galactic 7-ray diffuse emission is of 
utmost importance as this component needs to be sub- 
str acted to identify other extended astrophysical sources, 
such as the extra-galactic background or a possible dark 
matter component, and also because it is an interest- 
ing probe of the cosmic-ray population not only locally, 
but all over the Galaxy. In Delahaye et al. [9], we have 
shown how the morphology of the main component of the 
7-ray diffuse emission (namely the one due to 7r° decay) 
is affected by the various uncertainties due to cosmic-ray 
propagation, gas distribution etc. Here we propose to 
look more closely into the specific features of the spec- 
trum of this emission, especially in the light of the recent 
experimental data [1]. 

The local 7-ray emissivity from tt° decay has been 
thoroughly studied by the Fermi collaboration [1] which 
found that the data are consistent with a high nuclear en- 
hancement factor cm of 1.84 as found by Mori [17] when 
considering the cross-section given by Kamae et al. [13]. 
However we have already seen (see ref [9]) that these data 
seem to be in tension with the more recent cross-sections 
by Huang et al. [12]. 

In this letter we investigate explanations of this dis- 
crepency as well as the various uncertainties affecting the 
estimate of the 7-ray emissivity. We quantify the impact 
of the primary cosmic-ray spectra, of their propagation, 
of the metallicity of the Interstellar medium (ISM) and 
of the production cross- sect ions. For this, we follow the 
method developed in Delahaye et al. [9]. As a reference 
case, we will consider the primary fluxes by Shikaze et al. 
[24], the med propagation parameters, a helium to hydro- 
gen ratio of 1/9, and the cross-sections by Kamae et al. 
[13] with the nuclear weights from Norbury & Townsend 
[18]. After having defined the various quantities at stake, 



310 



Delahaye et al. 2011 
t 1 1 — i — i — I I I I 



> 
ID 



X 



5 10 




10 ] 



10° 

Energy [GeV] 



10 1 



FIG. 1. Effective emissivity as a function of energy. Plain 
lines correspond to Kamae et al. [13] cross-sections and dashed 
lines to Huang et al. [12] ones. Red and orange lines are 
respectively for med and min propagation parameters. The 
green line corresponds to the Kamae et al. [13] cross-sections 
used with a nuclear enhancement factor of 1.84 as suggested 
in the Fermi paper. The blue and purple lines correspond 
respectively to a minimum and a maximum case which are 
detailed in the text. Gray points are the recent Fermi data [1] 



we vary each parameter one by one to estimate its impact 
on the 7-ray emissivity and finally conclude. 
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DEFINITION OF THE QUANTITIES OF 
INTEREST 

The hadronic 7-ray emissivity of an interstellar nucleus 
A impacted by the cosmic-ray species b is given by the 
convolution: 



ei(x,E)-- 



da 
dE 



(&[T]+A^ 7 [£])x$ b (x,T) 



where 3>b(x, T) is the flux of cosmic ray species b. Be- 
cause the 7-ray flux is correlated to the hydrogen num- 
ber density nn, it is more convenient to use the effective 
emissivity per hydrogen atom defined as. 

£ eS (*,E)=J2— {£ P A (x,E)+£Z( X ,E)}. 
a nH 

The Fermi collaboration has selected two intermedi- 
ate latitude regions in which the hydrogen is atomic and 
mostly local. The point sources and the inverse Compton 
component have been substracted substracted from the 
maps. The residual flux has been correlated to the HI 
column density, yielding an effective emissivity £ e ^{E) 
which could be considered a priori equal to the solar 
value £q(E) = £ e ff (x , £"). How local is the Fermi mea- 
surement is at the center of our analysis. What Fermi 
actually measures is in fact quite close to the average 
of the 7-ray emissivity £ e ^ over the lines of sight of the 
pixels of the maps. What Fermi actually measures is in 
fact quite close to the average of the 7-ray emissivity £ e ^ 
over the line of sight. That is why we have computed 
the emissivity inside each pixel by weighing it with the 
HI spatial density and averaged the results over the sky 
region investigated in the Fermi analysis. 

A quantity commonly used in the literature is the nu- 
clear enhancement factor which is defined as the ratio 
= £ e ff(x, £ , )/f^( x 5 E)- The variations of cm with po- 
sition x and energy E have been so far disregarded and 
the nuclear enhancement factor has been mainly intro- 
duced as a constant by which proton-hydrogen interac- 
tions have to be renormalized in order to yield the total 
7-ray flux. Even at the Sun, it depends on the energy 
since it can be expressed as 



e M (E) 



n A f 



w(l,A) + w(4,A) 



*p(©,£) 



where, w(Ai,A2) are the nuclear weights by which the 
proton plus proton cross-section is multiplied to get the 
one of heavier elements collision. Untill recently, it was 
thought that the ratio of a to proton cosmic ray fluxes 
is constant with energy. However, recent results from 
[2, 21, 26] indicate that this ratio is increasing at high 
rigidity. The exact explaination of this increase is still 
under investigation (see for instance [4]) but has quite an 
impact on all secondary cosmic ray [11, 14] fluxes. As one 
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FIG. 2. Nuclear enhancement factor ploted as a function 
of the energy in the specific case of Kamae et al. [13] cross- 
sections with nuclear weights from Norbury & Townsend [18]. 
Other models give similar results. The colored curves corre- 
spond to primary cosmic ray fluxes measured by BESS [24], 
CREAM [26], ATIC [21] and Pamela [2]. 



can see from Fig. 2, cm may vary by more than 30% from 
1 GeV to 1 TeV making this quantity absolutely useless 
especially in the case where the fit to the PAMELA data 
proposed in the appendix are used. 



PROPAGATION 

As one can see from Fig. 3, in the direction I G 
[-160°; -100°] ,6 e [22°; 60°], which is the direction of 
interest of the Fermi study, the gas is mainly within 1 kpc 
from the Sun. However, depending on the propagation 
parameters chosen, the cosmic ray flux may exhibit a 
gradiant within this 1 kpc. This translates in the fact 
that £ e ff can differ from £@ by up to 10% especially at 
energies higher that 10 GeV (see Fig. 4). 

In the present case, we have taken into account the 
cosmic-ray propagation within the framework developed 
in Delahaye et al. [9] making use of the propagation pa- 
rameters which give a good agreement with all cosmic- 
ray flux measurements [16, 22]. Within this consistent 
framework three propagation parameter sets have been 
singled out [10], they give a relatively good estimate of 
a median and two extreme cases labeled min, med and 
max in Figs. 1,3 and 4. The med set is always used in this 
work unless specifically written otherwise. 

Note however that the gradient of cosmic-ray in the 
direction of study not only depends on the propagation 
parameters but also on the cosmic-ray source profile. As 
one can see in Fig. 4 if we consider the less steep source 



3 



r [kpc] Delahaye et al. 2011 

8.6 8.8 9 9.2 9.4 9.6 9.8 10 



Delahaye et at. 2011 



Gas distribution and cosmic ray gradient 
_{l,b) = {[-160 o :-100 o L[22 o :60°]) 




ur 4 • 

3 



io 5 



3= 



0.95 - 




s [kpc] 



10 1 
Energy [GeV] 



FIG. 3. Gas density (left ordinate) and cosmic ray gradient 
(right ordinate) as a function of distance from the Sun (lower 
abscissa) or from the galactic center (upper abscissa). 
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FIG. 4. Effective emissivity divided by solar emissivity as a 
function of energy. Plain lines correspond to source distribu- 
tions by Lorimer [15] and dashed lines to the one from Paczyn- 
ski [20]. Red, orange, and blue lines respectively correspond 
to med, min, and max propagation parameters sets. 



FIG. 5. Effective emissivity divided by solar emissivity as a 
function of energy. Various values of the helium to hydrogen 
ratios for the interstellar gas. In red the value considered as a 
reference in this paper, in pink a value consistent with obser- 
vation of the local region and in dark green a value consistent 
with observations of the anti-center. 



abundance of helium with respect to hydrogen is not pre- 
cisely known. The cosmological value is about 0.077 He 
atom per H [19] however, due to stellar evolution, the 
metallicity has evolved in the Milky Way. Because the 
stars do not have the same mass and age everywhere 
in the Galaxy, one should expect to have a helium gra- 
dient decreasing from the Galactic Center towards the 
outside regions. This is indeed supported by various ob- 
servations [5, 8] and theoretical models [6]. On average 
one can consider that nn e /^H varies from ~0.111 at the 
Galactic Center to ~0.087 in the outer region with a lo- 
cal value of ~0.097. These three values are illustrated in 
Fig. 5. The reason why we have considered a high metal- 
licity as a reference model is only because it was the value 
taken by Huang et al. [12], allowing easier comparisons. 

The p+He process and the a+He one amount respec- 
tively for ~20% and ~5% of the total signal hence vary- 
ing nHe/^H by 20% translates in a variation of the total 
signal of 5%, as one can see from Fig. 5. 



profile by Paczynski [20] rather than the one by Lorimer 
[15], the effect of propagation is quite reduced. 

METALLICITY 

The chemical composition of the interstellar medium, 
and what is more interesting in our case, the relative 



CROSS-SECTIONS 

As a reference case (hereafter called KNT) we con- 
sidered the cross-sections by Kamae et al. [13] (based 
on Pythia 6 [25]) with the nuclear weights from Nor- 
bury & Townsend [18]. Huang et al. [12] (but also Mori 
[17], which lead to the value of cm = 1.86) is based on 
DPMJET3 [23]. Both these Monte Carlo generators were 
not tuned with TeVatron and LHC data at the time of 
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FIG. 6. Effective emissivity divided by the reference so- 
lar emissivity (KNT) changing cross-sections and nuclear 
weights. 



7-ray production by 7r° decay is far from negligible. 

The statistical uncertainties of the Fermi results vary 
from ~4% to ~30% which is of the same order than all 
the effects we have discussed here. Moreover, the Fermi 
data stop at 10 GeV whereas most models sensibly dif- 
fer from each other at higher energies. It will hence be 
of great interest to have more data in the coming years. 
However, improving the experimental precision calls the- 
oreticians to produce as much an effort to reduce the 
theoretical uncertainties we have pointed here. 
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supported by the Spanish MICINNs Consolider-Ingenio 
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grant FPA2009-08958, the Community of Madrid under 
grant HEPHACOS S2009/ESP-1473, and the European 
Union under the Marie Curie-ITN program PITN-GA- 
2009-237920. 



these works hence it is probable that one may expect 
some change in the future. Though Cholis et al. [7] have 
shown that taking into account more recent versions of 
Pythia does not change much the results. 

However, considering what is available for the moment, 
as one can see in Fig. 6, the different models of nu- 
clear weights and cross-sections have a large impact on 
the estimate of the 7-ray emissivity. In the case of the 
production cross-sections of Huang et al. [12], the differ- 
ence reaches 50%. Note that the tables available concern 
cross-sections of cosmic-ray protons and a which interact 
with a mixture of hydrogen and helium but also carbon, 
nitrogen and oxygen. It is hence not possible to change 
these proportion and it is hard to tell whether the differ- 
ences come from the interstellar medium metallicity or 
from the Monte- Carlo. 

The sharp features of the green line of Fig. 6 are due 
to the resonnaces that Huang et al. [12] have added to 
their cross-sections. 



CONCLUSION 

The various effects we have pointed out in this letter: 
primary cosmic-ray fluxes, propagation, and metallicity 
which lead to uncertainties of order 45%, 10% and 5% re- 
spectively are illustrated in Fig. 1 in the KNT case by the 
blue and purple lines which show the two extreme cases. 
On top of that, we have also considered the uncertain- 
ties due to 7-ray production cross-sections and nuclear 
weights which vary from 10% to 50%, depending on the 
energy. Hence the total theoretical uncertainty on the 
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Fit to the Pamela data 

In order to do the computations we mainly used the 
fit performed by [24] over the BESS data but also the 



fit done by Lavalle [14] over data from ATIC [21] and 
CREAM [3] balloons. Moreover, we propose here a fit to 
the recent PAMELA data [2] : after demodulating in the 
force field approximation with a Fisk potential of 500MV: 

# P CT) = 35.310- (l - exp (- (^f 9 )) (to^v)" 2 ' 5 
, _r_x-o.5 (l T (1 _r_ r o.2i 

A I" 1 - ^ WGeVJ \ L T 300GeVJ V 1 ^ 5TeV/ 

and 

* a {T) =1.5 x 10- (sogv)" 2 - 7 (1 + db)" 13 

where R stands for the rigidity of the particle and both 
fluxes are expressed in ((GeV/n).s.sr.cm 2 ) -1 . 



